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SUMMARY 


The object of this report is to provide an estimate, based on 
theoretical calculations, of the forces induced on a wing that is fly- 
ing at a constant forward speed and suddenly enters a vertical gust. 

The calculations illustrate the effects of Mach number (from 0 to 2) and 
aspect ratio (2 to »), and solutions are given by means of which the 
response to gusts having arbitrary distributions of velocity can be 
calculated. The effects of pitching and wing bending sore neglected and 
only wings of rectangular plan form are considered. Specific results 
are presented for sharp-edged and triangular gusts and various wing-air 
density ratios. 


INTRODUCTION 


Studies of the gust-response problem for restrained wings (wings 
of infinite mass) entering sharp-edged gusts at supersonic speeds are 
already well advanced. Miles, Strang, Biot, and Heaslet and Lomax 
(refs. 1, 2, 3, and 4) presented solutions to such problems for two- 
dimensional wings; Miles and Goodman (refs. 5 and 6) presented solutions 
for rectangular wings having tip Mach cones that do not intersect the 
opposite edge. Miles and Strang (refs. 7 and 8) gave results for a tri- 
angular wing with supersonic edges. Theoretical studies restricted to 
incompressible flow fields contain the classical solutions due to Wagner 
(ref. 9) , Kussner (ref. 10), and von Karman and Sears (ref. 11), the 
former containing the solution for the indicial lift on a two-dimensional 
sinking wing and the latter two containing the solution for the lift on 
a restrained two-dimensional wing entering a sharp-edged gust. The 
extension of these studies to include the gust response for wings of 
finite aspect ratio has been carried out by Jones (ref. 12) . Later, 
further extensions to include the effects of gust shape as well as aspect 
ratio have been made by Zbrozek (ref. 13) and Bisplinghoff , Isakson, 
and O'Brien (ref. 14). 



2 


NACA TN 2925 


The purpose of the present report is twofold: first, to present 
solutions for a two-dimensional restrained wing entering a sharp-edged 
gust at sonic and subsonic Mach numbers (specifically, Mach numbers 
equal to 1.0, 0.8, and 0.5)j and second, to use these results together 
with those mentioned above to estimate the effect of wing aspect ratio 
and airplane mass on the lift response for airplanes flying at various 
speeds through the Mach number range from 0 to 2 and penetrating both 
triangular and sharp-edged gusts. 

A list of symbols is given in the appendix. 


ANALYSIS 

The Equation of Motion 


If induced pitching moments are neglected, the motion of a rigid 
wing disturbed from its equilibrium position by arbitrary external 
lifting forces is governed by Newton's second law. Thus, if w is the 
vertical velocity of the wing and m its mass, we can write 


(1) 

where the forces to be summed are the aerodynamic ones due to the gust 
velocity and the motion of the wing from its position of equilibrium. 

First, consider the force that results from a small vertical motion 
of the wing. Suppose the wing has been flying in steady level flight at 
a constant speed U Q up to a time t' =0. Fix an xyz coordinate 
system in space (z positive upward) such that at , t' =0 the y axis 
lies along, and the origin on, the wing leading edge and, further, such 
that the wing is moving in the negative x direction. For t’>0 the 
wing moves away from these coordinates, continuing forward at the con- 
stant speed U 0 along the negative x axis , and now also moving down- 
ward at a constant rate -w = U Q a. The transient lifting force on the 
wing induced by such a maneuver shall be referred to as the indicial 
lift (positive upward) and designated in coefficient form by the sym- 
bols c Z a or f° r section- or total-lift values, respectively. ' 

Given, the indicial lift coefficient, one can show by using the 
principles of superposition that the lift due to an arbitrary variation 
of angle of attack caused by the vertical velocity of the wing can be 
determined from the relation 

L(t ' ) =_ UodF f c iu (t ’-VMV) dt i 

0 J o 


dw V 

m = ) 

dt' 


forces 


(2a) 
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Next consider the force that is 
induced on a wing penetrating a 
sharp-edged gust having a uniform 
upward velocity w D . If the wing is 
restrained so that it can move 
neither upward nor downward (corre- 
sponding in flight to the limiting 
case of infinite wing mass), the 
section or total lift coefficients 
induced by a unit value of w 0 /U 0 
will be designated ct or Cl , 

respectively. Sketch (a) illus- 
trates the differences between the 
boundary conditions for, and the 

initial variations of ci 

u cl 

and ci for a two-dimensional wing 
traveling at a subsonic speed. 


Given a value of one 

obtains the lift on a restrained 
wing flying into a gust having an 
arbitrary vertical velocity dis- 
tribution, w a (t’), by the relation 


Section 

leading 

edge 

trace 



Section 

trailing 

edge 

trace 


w - - U 0 a 





Sketch (a) 


r' < 2b > 

U Q dt' & 

Substituting equations (2a) and (2b) into equation (l), one finds 
the expression for the vertical motion of an unrestrained wing flying 
into a gust; thus 


m 


r l Ct (t'-t ')w(t')dt . 

dt' U Q dt' J ^ 11 1 

o 


(3a) 


Since w a (t') is assumed to be given, equation (3a) is an integral equa^ 
tion - in terms of w(t') - of the second kind with a variable upper 
limit. It is convenient first to study equation (3a) when the gust is 
a step function (sharp-edged gust) • For this case w a ( t ' ) becomes a 
constant w Q , say, and equation (3a) reduces to 


m 


dw 

dt' 


Ct (t'-t 1 ')w(t 1 ')dt ' 

U n dt' J ^ 

u o 


+ 5&p, (f) 


u, 


(3b) 
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The solution to equation (3b) can be used to find the induced force 
on an unrestrained wing entering a gust of arbitrary structure. Methods 
lor solving the integral equation and applying its solution will be 
developed in the subsequent sections. 


Indie ial Lift on a Sinking Wing 


The analysis involved in calculating the indicial lift force on the 
wing is based on the assumptions that underlie linearized, thin-airfoil 
wing theory in general. Mathemetically, these assumptions imply that the 
governing partial differential equation of the flow field is the three- 
imensional wave equation. In terms of the velocity potential cp and 
for an axial system fixed relative to the still air at infinity, the 
wave equation can be written 


'Pxx + 'Pyy + <P Z z ~ <Ptt 


(4) 


where t is the product of the speed of sound a Q and the time t*. 

For a wing moving in the z = 0 plane, the boundary conditions are that 
T is continuous everywhere except across the wing and its vortex wake, 
^z) z=0 is a constant over the region bounded by the wing plan form at 
given time, and all velocities vanish outside the starting wave 

PTlVPl rvn»=k o . 


All values of Cj^ and c^ used herein have been presented in 

previousiy pubiished reports. The indicial lift on a sinking rectangular 
wing traveling at supersonic speeds has been presented by Miles (ref? 5). 
As the aspect ratio tends to infinity, this solution approaches that for the 
two-dimensional case given in references 1, 2, 3, and 4. At sonic and 
subsonic speeds results, for a two-dimensional wing are available for 
Mach numbers equal to 1.0, 0.8, and 0.5 (ref. 15) and for the incom- 
pressibie case (ref.l 6 ). Finally, the indicial lift on sinking wings 
finite aspect ratio in incompressible flow is presented in refer- 
ence 12. Curves showing the effect of Mach number on the two-dimensional 
values are. presented in figure l(a), and the effect of aspect ratio at 
supersonic speeds is indicated in figure 1(b). Tabular values for the 
two-dimensional wing flying at Mach numbers equal to 0.5 and 0.8 are 
given in tables I and II, respectively. 


Response of a Restrained Wing to a Sharp-Edged Gust 


Load distribution.- The lift induced on a restrained wing penetrating 
a sharp-edged gust can also be determined by solving equation (4) subject 
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to the proper boundary conditions. 1 For a wing moving in the z = 0 
plane these conditions are similar to those given for the indicial lift 
on a sinking wing in that all velocities vanish outside the starting 
wave envelope and 9 is continuous everywhere except across the wing 
and its vortex wake, but differ from the indicial case in that (<p z ) z _q 
is a constant only over the portion of the wing plan form that has pen- 
etrated the gust, being zero over the remaining portion (see sketch (a)). 
This problem has been solved for a rectangular wing traveling at super- 
sonic speeds by Miles (ref. 5) and, again, as the aspect ratio tends to 
infinity, this solution approaches that for the two-dimensional case 
given in references 1, 2, 3> and 4. Two-dimensional wings flying at the 
speed of sound and two- and three-dimensional wings flying in an incom- 
pressible medium have also been considered (refs. 4, 16, and 12, 
respectively) . 

The problem of finding the two-dimensional gust response at subsonic 
speeds can be solved by the same method that was used in reference 15 to 
find the two-dimensional subsonic indicial response. For these cases 
equation (4) reduces to 

^xx + ^zz = ^tt (5) 


and the boundary conditions for a section in the xt plane are indicated 
in sketch (b) 


~ c ~j 

Trace of foremost sound /% 

wave ✓ M — t 2 * 

/-& §_J^ 

Loading falls to zero / /pi- / 
at trailing edge — 7^ — 

Vortex wake, A 6 Jl|§ 

independent of t Jlilff 

and given by value lltffl 

at trailing edge jftllfl-ifl 



Sketch (b) 


It is interesting to note that the gust lift function c j can be 

D 

related to the indicial response following a step variation of angle 
of attack under quite general "conditions by the reciprocal theorems 
given in reference 17 . 
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The solutions obtained for the load coefficient over regions 1 and 2 
shown in the sketch can be written (details of the analysis are omitted) 

For region 1 


AP 

qg 



8 /Mo(t-x) 

ir(l+Mo) J x+Mot 


For region 2 


AP 

qg 


■jt(l+Mo) I + Jt '/ M o( t “ x )( c - x - M o t ) 


2K 


' EF'(t)+KE'(t)-KF , (t) 

7 (t 2 -x 2 ) (1 -Mq 2 ) y ( x+Mot ) ( c-x-Mot) 


(6a) 


(6b) 


The symbols E,K,E'(ty) and F'(^) are elliptic integrals defined in the 
appendix, their modulus k being given by 


k = 


I- 


(t+x)(l+Mo)-2c 
( t+x) (1 +Mq) 


and their argument ^ by 


(7) 


. . /x+Mot 

V = arc sm J 


(8) 


Equations (6a) and (6b) give the loading over the complete wing 
section for values of t less than or equal to 2c/(1-Mq 2 ). Since 
M 0 t/c equals U 0 t'/c, the number of chord lengths traveled, these 
equations represent the exact linearized solution for the section load 
distribution during the time required for the w ing to travel 
2Mo/(l-Mo 2 ) chord lengths after reaching the gust front. Hence, for a 
Mach number equal to 0.8, equations (6) establish the gust response dur- 
ing the. first 4.44 chord lengths of travel after the gust penetration. 
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Sketch (c) shows the variation of 
gust loading Ap/qg throughout 
this interval and also, for com- 
parative purposes, the indie ial 
load variation Ap/qa for the first 
k chord lengths of travel. The 
dashed curves in the sketch repre- 
sent the final steady- state load 
distribution adjusted so as to give 
the same total lift as the exact 
solutions for the gust and indicial 
cases at U 0 t'/c equal to 4.44 d 
and 4, respectively. Thus, to the - 
degree of accuracy indicated in the 
sketch, the gust and indicial load- 
ings at Mq = 0.8 can be approxi- 
mated for larger values of U 0 t’ /c 
by the expressions 


c Z g (Uo 

c Z g 0 


it'/c) 

—j 



Percent chord 


2Mq U Q t« 
1-Mo 2 c 

and 


(9a) 



The variation of c 7 and c, 

, g L a, 

for values of Upt'/c greater 
than 4 will be discussed presently. 

For a Mach number equal to' 0.5, 
equations (6) are .sufficient to 
establish the gust response for 
J?nly f ir st_1^33 cho rd l engths _ 
traveled. Further calculations 
were carried out and the exact 



Sketch (c) 
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loading was established for both 
the gust and indie ial cases for 
values of U 0 t'/c less than or 
equal to 2.33* These calculations 
were for the most part numerical 
and no simple closed expressions 
such as those presented in equa- 
tions (6) were obtained. Sketch (d) 
contains the results. Again, the 
dashed curves represent the final 
adjusted steady-state load distri- 
bution indicating that the gust and 
indicial loadings for Mq = 0.5 can 
also be approximated for larger 
values of U 0 t'/c by the equa- 
tions ( 9 ). 


Percent chord 


too 


kc-H* 

* * 


Section lift .- When integrated 
across the chord at a fixed time, 
the loadings shown in sketches (c) 
and (d) give the variation of the 
lifting force on the wing section 
during the early portion of the 
response. In the interval 
0 ^U Q t' /c <Mq/(1+Mq) equation (6a) 
integrates to give 


c l 

V U o 


% = 



( 10 ) 



In the interval 

M 0 /( 1+M q ) < U Q t ' /c < 2M 0 /( 1-M 0 2 ) 

the expression for the loading is 
too complicated to integrate analyt- 
ically and the section lift was 
calculated by numerical methods. 

The results , together with those 
for c? (taken from ref. 15 ), are 
given in sketch (e). Since, as time 
goes on, C 2 g must approach cj a , the 
curve for the gust response was simply 
faired into the curve for the indicial 
response in the manner shown by the 


Sketch (d) 
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dashed lines. Finally, for values of U 0 t’/c greater than 10, the 
following equations , taken from reference 15, can be used: 

For Mq = 0.5> 

~ -2« r x _ 4 1 44.218 


c lg ~ c la 


JT- 


Mq" 


3 5+2(U 0 t'/c) [5+2(U 0 t'/c)]‘ 


} 


(11a) 


F or Mq — 0 • 8 y 


„ 2* I 1 1.73b TO. 83 1 

c lg ~ c la~ /n _ M 2 \ n+1.2^ U_t* /c) [ll+1.25(u 0 t' /c) ] 2 J 




1.736 

ll+1.25(U 0 t'/c) [11+1. 


(lib) 


The final curves (determined from the previous analysis and afore- 
mentioned references) for c lg> 4he section lift coefficient developed 
on a restrained wing entering a sharp-edged gust, are shown in fig- 
ure 2(a) for Mach numbers equal to 0, 0.5, 0.8, 1.0, 1.2, 1.41, and 2.0. 
Tabular values are given in tables I and II for Mach numbers equal 
to 0.5 and 0.8. The effect of aspect ratio at a Mach number equal to 1. 
is shown in figure 2 (b) . 


Response of an Unrestrained Wing to a Sharp-Edged Gust 


Given the indicial lift response C]^ and the response for a 
restrained wing penetrating a sharp-edged gust C Lg , one can use 


UJ*. 
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equation (3) to find the motion of an unrestrained wing entering a sharp 
edged gust having a constant upward velocity w 0 . As in reference 4, 
the lift on the unrestrained wing can be related to an infinite series 
of integrals involving and Ci g . First, set 


U 0 t« 

T = , 

C 


U 0 ti’ 


w 



2m 

P 0 cS 


so that, by integrating equation ( 3 ) with respect to t' , one finds 


W " 1 7 f C Lg( T l) dT l + zf C L a ( T - T l) w ( T i)dT 1 = 

o J o 


0 


Then use the relation 


( 12 ) 


CL 

w o/ u o 


= (dw/dr) (n/w Q ) 


and iterate equation (12) using Liouville’s method of successive sub- 
stitions. (See ref. 18.) This yields 


C L 

w o/C c 


1 r T 

= C Lg< T ) -- J c L a ( T -'l)CL g (T 1 )dT 1 +• 



T -T x )dT 1 



T i- T 2 ) c L g ( T 2 )<iT 2 - ... 


( 13 ) 
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Equation (13) converges uniformly 2 
for all T. By means of it, Cp/(w 0 /u 0 ) 
and cj/(w 0 /U 0 ), the total and section lift 
coefficients induced on unrestrained wings 
entering a sharp-edged gust, have been calcu- 
lated and the results are shown in figure 3* 
Chart (a) indicates the range of Mach numbers 
and aspect ratios for which calculations were 
made, the numbers in the chart referring to 
the individual figures in which the results 
are presented. It should be noted that 
results are given for a wing flying at Mq= 1 
and having a finite aspect ratio. Such 
cases can be calculated from the indie ial 
and restrained gust responses presented in 
reference 5* These responses are still 
valid at Mq= 1 for values of the time 
variable up to that for which the wave enve- 
lope induced by one side of the wing crosses 
the opposite side. An aspect ratio 5 wing 
flying at the speed of sound travels 13 chord 
lengths during this time interval, and this 
is sufficient to establish the significant 
part of the response to a sharp-edged gust 
for n < 300. 


L02\ 


t.20\ 


The chart also shows that the gust it. 4/ 31 I oj I or 1 

response for the unrestrained wing was I • 1 I I 

calculated at Mq= 0 for an infinite- Chart (a) 

aspect-ratio wing (for comparative purposes) 

but not for finite-aspect-ratio wings. The gust response on both 
infinite- and finite-aspect-ratio wings in incompressible flow have been 
studied extensively by means of operational methods in references 12, 13 , 
and 14. Where comparisons can be made, the results obtained in this 
report using equation ( 13 ) agree well with those given in the references 
mentioned. 


3i 3e 


2 The statement made in reference 4 on the convergence of this series is 
unnecessarily restrictive. Since the greatest values of C]^ 
and CLg are Cp^®) and Cp f ®) , the general term of the series does 
not exceed ^ 


[Clo/*)] 0 c L g (°°) / dT i y 1 dT 2 J 2 


that is, does not exceed 


— [Ct (“)] n Ct (») — T n 

u n La/ ;J L g v ' nl 


and by the ratio test the series converge uniformly. 
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Variation of gust 
intensity 

Path of point on 
wing 

Variation of lift 
on wing 



Response of Unrestrained Wing to 
Arbitrary Gust 


The function Cl/(w 0 /U 0 ) presented in 
the previous section can be thought of as 
the indicial gust response for lift on an 
unrestrained wing. In this sense it is 
apparent that the lift on a wing penetrat- 
ing a gust in which w is a function of 
the chord lengths traveled can be calcu- 
lated by superposition and is represented 
by the integral 


= — f T V ( T - T l) dT 

" iT -J 0 V“o u o 1 


(Ih) 





Sketch (f) 


By means of equation (lk), the lift 
induced on a wing moving at the constant 
speed U Q and entering a gust, the ver- 
tical velocity of which starts at zero 
and increases linearly with distance of pen- 
etration, is simply the integral of 
Cl/(w 0 /U 0 ). Thus, representing the sec- 
tion lift coefficient developed by a wing 
entering a gust with a unit gradient by 
the symbol Cl s , we can write 


c l s < t > 


r ^ dTl 

J "o/Uo 


(15) 


If the wing flies into a gust with a 
triangular- shaped distribution of w, 
having its maximum intensity w t a dis- 
tance h chord lengths from the front, it 
follows at once that the resulting lift 
response Cl/(w-(-/U 0 ) is given by 


cl 

V U o 


I [C Ls (t)]j 0<T<h 
( i t C Lg ( t ) - 2C Lg ( t - h) ] ; h<T<2h > 

i [ C Lg ( r ) - 2C Lg ( t - h) + 

C L (r-2h)]; 2h<T 

L 3 

(16) 
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Examples of the various gust shapes and the responses in 
vertical motion of wings penetrating them are shown in sketch 


lift and 
(f). 


MAXIMUM LIFT DUE TO GUST PENETRATION 
Sharp-Edged Gust 


Consider the maximum increase in lift caused by the entry of the 
wing into a sharp-edged gust. This increment is given for the range of 
Mach numbers, aspect ratios, and wing-air density ratios shown in 
chart (a) by the maximum values of Cl/(w 0 /U 0 ) and C 2 /(w 0 /U 0 ) on the 
curves shown in figure 3* 

First let us consider wings of infinite aspect ratio. For such 
wings the variation of the maximum gust-induced lift coefficient with 
Mach n umb er is shown in figure 4, and a cross plot in which Mq instead 
of p is held constant is presented in sketch (g). The values 



Sketch (g) 
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V 

00 

300 

200 

0 

0 

0 

0 

0.8 

67 

39 

37 

1.0 

00 

\/03 

88 

1.2 

-4 

to 

12 


Percent increase in 
relative to 

max 

its value at Mo = 0 
Chart (b) 


for [x 00 are the steady-state values given 
by the simple equations 

WI - Mo 2 ; Mq< 1 

V <s/Mq 2 - 1; Mo>l 

The difference between the lift increment on 
a restrained wing and that on one with a 
finite value of |i is seen to be most pro- 
nounced at the high subsonic Mach numbers. 
Notice, for example, that the percentage • 

increase in [c^/(w 0 /U Q )] found by 

max 

increasing Mq from 0 to 0.8 is 67 for q=» 

( Prandtl-Glauret rule) but only 37 for p=200. 
Chart (b) indicates the relative increase in 
t c z/( w o/ u o)J inax caused by compressibility for 

three different values of the wing-air density 
ratio. J 



Consider next the effect of aspect ratio on the maximum lift 
increment induced on a rectangular wing penetrating a sharp-edged gust. 

+h ! 1 •?+" °* thl ? ir * crenfint is again given by the steady-state value of 
the lift-curve slope and is presented for A = 00, 5 , and 2, in figure 5 

liftin Stead r StSt J Values are taken from the numerous studies made of 5 ’ 
g sur aces traveling at subsonic and supersonic speeds. On the 

supersonic side, for cases in which 
e.g., ref. 19 ) 


a7V^7 — 1 , the equation (see. 


' 'Cl \ = 4 

' w o/ U o •'max J Mq 2 - 1 



applies and, for cases in which > 0 , the curves in 

reference 20 were u sed. On the subsonic side the portions of the curves 
m the range . §<kj\ - M Q 2 <2 were again taken from reference 20. 

The results in references 19 and 20 are sufficient to cover the entire 
Mach number range for the A = 2 wing. For the A = 5 wing the values 
on the subsonic side outside the range 0<kj\ - Mq 2 <2 were taken 
from a curve 3 that was compiled from a large number of solutions for 
lifting surfaces traveling at subsonic speeds. 


The curve was taken from an article prepared by Robert T. Jones and 

Doris Cohen for the forthcoming series on High-Speed Aerodynamics 
and Jet Propulsion, Princeton University Press. 
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The values of [C L /(w 0 /U 0 )] mfly for rectangular wings traveling at 

supersonic speeds, given in figure 3, and the incompressible-flow 
solutions, given in references 12, 13* and l4, were used to prepare the 
the curves in figure 6. The dashed lines between the Mach numbers of 0 
and that for which A*/Mq 2 - 1 = 1 are interpolated, the two-dimensional 
results presented in figure 4 lending credence to the validity of the 
interpolation. Sketch (h) presents the aspect-ratio effect on 
[C L /(v 0 /U o )l at Mq = 1. 



It should be noted that in 
the vicinity of Mq = 1, the 
curves for which A = p = 00 
( figs . 4 and 5) > and probably 
also those for which A = 5* 

= 00 (figs. 5 and 6(a)), are 
not valid representations of 
the gust-induced lift on actual 
wings flying at these speeds, 
although they do represent 
solutions to equation (4) con- 
sistent with the boundary con- 
ditions previously discussed. 

For the two cases, mentioned, 
the assumptions on which equa- 
tion (J+) is based axe violated. 

These assumptions are more 
closely approached, however, 
as the wing-air density ratio 
and the aspect ratio decrease. 

Hence, for lower values of A 
and p the solutions given herein for wings traveling in the transonic 
speed range have justification on a physical as weil as a mathematical 
basis. 


Sketch (h) 


Triangular Gust 


The maximum increase in lift on a two-dimensional wing passing 
through a triangular gust having its maximum intensity 12 chord lengths 
from its front is shown in figure 7- For the lower values of p the 

variation of [c 7 /(w + /u_)] with Mach number is similar to that cal- 

L ' 0 max 

culated for the sharp-edged gust and shown in figure 4. As p 
increases, however, a comparison of the results shown in these two 
figures indicates the - importance of the assumed gust shape in estimating 
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\type 
M \ v 

M e =0.8 

M o =/.0 

Sharp 

edge 

Tri- 

angle 

Sharp 

edge 

Tri- 

angle 

60 

1.24 

1.24 

1.44 

1.51 

100 

1.29 

1.29 

159 

1.58 

200 

1.37 

1.33 

1.88 

1.67 

300 

1.39 

1.34 

203 

1.70 

OO 

1.66 

1.36 

OO 

1.85 


Ratio of the value of (CL) max 
at Mq = 0.8 and 1.0 to its value 
at Mq = 0. 

Chart ( c) 


the maximum gust-induced lift. 
Chart (c) shows the difference 
in the compressibility effect 
obtained for the sharp-edged and 
trisingular (12 chord lengths to 
apex) gusts . 

Figure 8 presents the aspect- 
ratio effect on the maximum lift 
response for rectangular wings 
penetrating triangular-shaped 
gusts. The values at M 0 = 0 
were calculated from the results 
given in reference 12 and again 
the dashed lines represent an 
interpolation . 


CONCLUDING REMARKS 


Results are presented for the 
lift developed by a restrained 
two-dimensional wing flying at a 
Mach number equal to 0-5 or 0.8 
and penetrating a sharp-edged 
gust. Similar results are reviewed 
for Mach numbers equal to 0, 1.0, 
1.2, 1.41, and 2. 


A method is given whereby the lift can be estimated (neglecting the 
effects of airplane pitching and wing bending) for unrestrained rectan- 
gular wings in the aspect-ratio range 2 to <*>, flying in the Mach number 
range 0 to 2, and penetrating gusts of arbitrary structure. Specific 
results are given for sharp-edged and triangular-shaped gusts . 


In general, given variations in the wing aspect ratio, the wing- 
air density ratio, and the gust shape have their maximum effect on the 
gust lift when the wing is flying at a high subsonic speed. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Feb. 3, 1953. 
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APPENDIX 
LIST OF SYMBOLS 


A aspect ratio 

a Q speed of sound 

c chord length 

Cl wing lift coefficient, L/qS 

cj section lift coefficient, L/qc 

E complete elliptic integral of second kind with modulus k 

E'(i|0 incomplete elliptic integral of second kind with modulus k' 
and argument- i|r 

Y'(if) incomplete elliptic integral of first kind with modulus k’ 
and argument i|f 

g w Q /U 0 

h number of chord lengths from front to apex of triangular gust 

K complete elliptic integral of first kind with modulus k 

k modulus of elliptic integrals 

(See equation (7)0 

k» J 1 - k 2 

L lift on wing 

Mq Mach number at which wing is traveling 

m mass of wing 

Ap/q loading coefficient, pressure on lower wing surface minus 

pressure on upper wing surface divided by dynamic pressure 


dynamic pressure, i p 0 U Q 2 


S wing area 
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t* time 

t a 0 t* 

U Q wing velocity 

w vertical velocity of wing 

w a(t' ) velocity of arbitrary gust 
w o velocity of uniform, sharp-edged, gust 

w t maximum velocity of triangular gust 

x,y,z Cartesian coordinates fixed with reference to still air at 
infinity, z positive upward, y parallel to wing leading 
edge, negative x direction corresponding to direction of 
wing motion 

a wing angle of attack 

d wing-air density ratio, 2m/p 0 cS in analysis of complete wing, 

2m/Po c i n analysis of wing section 

P Q air density 

t chord lengths traveled by wing, U Q t' /c 

<p perturbation velocity potential 

^ argument of elliptic integrals 

(See equation (8).) 

Subscripts 

g response of restrained wing to unit, sharp-edged gust 

s response to gust with unit velocity gradient 

a indicial response on sinking wing 
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TABLE II.- VALUES OF c 7 AND c 7 FOR M n = 0.8 

£, a ‘'g 





U 0 t’/c 

yi-M o 2 =t a /2" 

y^Mo 2 c Zg/ 2 * 

0 

0.478 

0 

.1 

.466 

.044 

.2 

.454 

.085 

• 3 

.442 

.129 

.4 

.430 

.170 

•5 

.423 

.209 
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.426 

.234 
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.433 

.256 
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.442 

.276 

• 9 

.451 

.296 

1.0 

.461 

.315 

1-5 

.507 

.402 
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.546 

.465 
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.581 

.513 

3-0 

.610 

.551 

3.5 

.632 

• 584 

4.0 

.652 

.616 

4.5 ' 

.670 

.642 

5-0 

.687 

.663 

6.0 

.714 

• 700 

7.0 

.738 

-730 

8.0 

.760 

.758 

9.0 

• 779 

.780 

10.0 

.798 

.796 

00 

1 .000 

1.000 
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(a) Two-dimensional wing for several Mach numbers. 


Figure 1,- Variation of indicial lift response with chord lengths traveled. 
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Chord lengths trave!ed J U 0 t'/c 

(a) Mq = 0, A = » 

Figure 3*- Response of unrestrained wing to a uniform, sharp-edged gust. 
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Chord lengths traveled ,U 0 tyc 
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Q) 


Chord lengths t rave led, UJt'/c 
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Chord lengths traveled, U 0 t'/c 
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Chord lengths traveled \ Uj'/c 

(g) Mq = 1.00, A = 5 . 
Figure 3 .- Continued. 
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Chord lengths traveled, U 0 t'/c 

(h) Mo = 1.02, A = 5 
Figure 3 .- Continued. 
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Chord lengths traveled, U 0 t'/c 
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Mach number, M 0 

Figure 4.- Maximum increment, of lift induced on a two-dimensional wing 
entering a uniform' sharp-edged gust < 
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Mach number , M c 

Figure 5 .- Maximum increment of lift induced on restrained rectangular 
wing entering . a uniform sharp-edged gust . 




Figure 6.- Maximum increment of lift induced on a rectangular wing 
entering a uniform sharp-edged gust. 
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Figure 7;- Maximum increment of lift induced on a two-dimensional, 
wing entering a triangular gust having its apex 12 chord lengths 
from the front.. 
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NACA-Langley - 4 - 16-53 - 1000 


Figure 8.- Maximum increment of lift induced on a rectangular wing 
entering a triangular gust having its apex 12 chord lengths from 
the front . . 
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